Complex Energy Method in four Body Faddeev-Yakubovsky Equations 
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The Complex Energy Method [Prog. Theor. Phys. 109, 869L (2003)] is applied to the four body 
Faddeev-Yakubovsky equations in the four nucleon system. We obtain a well converged solution in 
all energy regions below and above the four nucleon break-up threshold. 
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Calculations for scattering systems in configuration 
space require boundary conditions which increase in com- 
plexity with growing particle numbers. These bound- 
ary conditions appear in the form of Green's functions 
in momentum space which carry singularities of increas- 
ing complexity. The Green's functions are expressed as 
Go = 1/(E + ie— Ho) where E and Hq are the total and 
kinetic energy, respectively, and the limit e — > has to 
be taken. In the two-body system there is one (relative) 
momentum variable p and Go has a pole in the complex 
p-plane. It is easy to handle it using the principal value 
prescription and (half) the residue theorem (PVR). In 
the three body case there arises already a difficulty in 
the form of so called moving singularities however, 
PVR is still applicable Q, or one can use the contour 
deformation |j, la, Ifl 13 (CD) technique. Summarizing 
these techniques, first one takes the limiting value e — > 
and next the equation is solved avoiding the integration 
path on the complex plane. This is illustrated in Fig. ^ 
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The situation is more complicated in the four body 
system. Employing a separable potential and a sepa- 
rable expansion technique for the three body, and [2+2] 
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FIG. 1: Illustration of integration paths for PVR (left) and 
CD (right). The cross in the left figure indicates a fixed pole. 
Moving singularities occur in the shaded area of the right 
figure. 

subamplitudes, the four body Faddeev-Yakubovsky (FY) 
equations Q for four identical particles can be expressed 
as 
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where the A4's are the four body amplitudes, a and 
indicate [3+1] and [2+2] configurations (see Fig. £ is 
an exchange term from [3+1] to [3+1] configurations of 
which the plus sign corresponds to the four bosons and 
minus sign to the fermions. ^'s are exchange terms from 
[2+2] to [3+1] where the subscripts are related to the two 
diagrams in Fig. 3. More details may be found in Ref. 0- 
Further H and Q are the three body and [2+2] propaga- 
tors and they have a similar nature as Green's function 



in the two body Lippmann-Schwinger (LS) equation. 

If one stays below the three-body break-up threshold 
the FY equations can be solved with PVRj_since only 
two-body singularities occur [l(J [U E E 111 E E] . 
Above the three-body threshold but still below the four- 
body threshold the FY equations have also been solved 
applying the CD techniques E E E E E E- 
There in the £ and J-'s terms occur two-body propagators 
whose nature is similar to the three body Green's func- 
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FIG. 2: The two partitions in four particle system, k, p's, 
and g's are standard Jacobi momenta. 
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FIG. 3: Diagrams for the £, Ti, T2, Ti., and Q ingredients 
to Eq. (1). The shaded parts with r indicate the two body 
propagators. 



tion in the Born term of the three body Alt-Grassberger- 
Sandhas 0] or Amado-Mitra-Faddeev-Lovelace (e.g. 
[2(j) equations. However, above the four-body thresh- 
old the four-body Green's function depends on all (rel- 
ative) momenta and the behavior of those singularities 
is quite complicated. Thus neither PVR nor CD tech- 
niques have been successfully extended at energies above 
the four body break-up threshold and we are not aware 
of a solution in this energy region. 

Recently the Complex Energy Method (CEM) has 
been revived and applied to the two- and three-nucleon 
system. The first step of CEM is to solve the equation 
with some finite e's. These calculations are easily carried 
out since there are no singularities on the real momentum 
axis. After obtaining solutions with various e's, the lim- 
iting value e — > is taken numerically with an analytical 
continuation method. 

Our aim is to generate solutions applying this method 
to the FY equations in all energy regions including en- 
ergies above the four body break-up threshold. We per- 
formed calculations in the J 77 = + and T = state for 
the four nucleon system. For this feasibility study the 
J OT = l/2 + state is included in the three body subsystem 
and the 1 So and 3 Si- 3 Di states in the two body sub- 
system. All allowed spins and angular momenta within 
this restriction are included ; in short, there are 14 chan- 
nels. The Coulomb force is neglected. The Yamaguchi 
potential [23, H3 is employed as the nucleon-nucleon in- 
teraction. The potential V has a separable form as 

V U '(k,k')=g e (k)Xg e/ (k'), (2) 

where the g£S are the two body form factors for the par- 



tial waves I as 

1 Ck 2 
9o(k) = 1 ^ TW , 92(h) = (fc2 + /?2)2 , (3) 

and the fc's are the initial and final momenta between 
the two nucleons, respectively. We adopt the parame- 
ters A, (3, and G given in Tabic |U We represent the 
three body and [2+2] subamplitudes by rank-4 separa- 
ble forms employing the Energy Dependent Pole Expan- 
sion [24| (EDPE) method. We take the nucleon mass as 
938.91897MeV which is the average of those for proton 
and neutron, and he = 197.327054 MeVfm. The inte- 
grations are cut off at 200 fm _1 for k, at 40 fm for the 
p's, and at 16 fm -1 for the q's (see Fig. [3J. 

The FY equations are solved at four energies: (I) 
1.5MeV above the 3N+N threshold, (II) lMeV above 
the 2N+2N threshold, (III) lMeV below the four body 
break-up threshold, and (IV) 12MeV above it (see Fig.0J. 
We define the ie term of the four body Green's function 
as is + where e and £ are real. Thus Go turns into 
Go = 1/ {E + ie + Q — Hq). Solutions of the FY equation 
satisfy uniqueness even at the limit for e —> 0, which is 
not the case for simple LS equation [2{|. Therefore the 
results by the analytical continuation do not depend on 
the choice of e's within the radius of convergence. Thus 
we empirically choose 0.5MeV as the minimum e value 
for the cases (I)- (III) and 0.75MeV for the case (IV) (see 
crosses in Fig.QJ, with attention only to a beter numer- 
ics. They are increased in steps of 0.125MeV. £ is chosen 
as and ±0.125MeV. 

We employ the point method [2(| as an analytical con- 
tinuation technique in CEM. It's convergence behavior is 
shown in Table ITU where the phase shift 8 and the inelas- 
ticity parameter rj is defined by S = 7yexp(2z<5). Here S 
is the S-matrix of elastic 3N+N scattering and is related 
to the on-shell amplitude of M aa in Eq. (JIJ (A4™ ) as 
S = 1 — 2i«;A / t°Q where re is the on-shell momentum. 

In case (I) r\ must be 1 due to unitarity and our result 
satisfies it within six digits. Also in cases (II) and (III) we 
reach a very high accuracy. In case (IV), we still obtain 
converged solutions within 4 digits. In the cases (I and 
II) our results agree very well with the solutions based 
on PVR. 

We showed that well converged solutions of the FY 
equations are obtained in all energy regions. In relation 
to the application of EDPE we confirmed that converged 
solutions are obtained in the cases (I-III) . In the case 
(IV), however, there is a report that EDPE is not appli- 
cable [2]]]. We also found that EDPE did not converge. 
Therefore, in this study we just kept the rank fixed by 



TABLE I: Parameters of the Yamaguchi potential. 
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A (MeV fm" 1 ) 


state 


P (fm- 1 ) C (fm ) 




-68.942626 


X So 


1.1300 


3 Si- 3 D! 


-74.506955 




1.2412 
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1.9476 -4.4950154 
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FIG. 4: Illustration of threshold energies for the 4N system. We choose E=0 at the four-body threshold. The various energies 
for the calculations are measure relative to the thresholds. The crosses indicate the complex energies where we solve the FY 
equations in the CEM. They are numbered by n for each choice of the energy region (I- IV). 



TABLE II: Phase shifts (in deg) and inelasticity parameters for 3N+N — > 3N+N elastic scattering. n max denotes the number 
of sample energies which are included in the point method. For instance, n max = 5 means that the solutions from n—1 to 5 
(see Fig. are included. The row PVR shows results from a direct solution of the FY equations using PVR. The agreement 
is perfect. 





(I) 


(II) 


(III) 


(IV) 


^max 


5 


V 


5 


n 


5 


V 


s 


n 


l 


53.31351 


1.018810 


14.62234 


0.675088 


-8.0617 


0.692391 


-62.093 


0.83875 


2 


46.21307 


0.973139 


10.47942 


0.853658 


-5.9428 


0.815698 


-61.965 


0.75946 


3 


44.27898 


0.989232 


12.38204 


0.948787 


-5.5022 


0.899150 


-61.620 


0.74499 


4 


44.27129 


1.000623 


12.38254 


0.948044 


-5.5101 


0.898666 


-61.676 


0.74570 


5 


44.34441 


0.999787 


12.38211 


0.948046 


-5.5095 


0.898656 


-61.682 


0.74589 


6 


44.34157 


0.999994 


12.38284 


0.948053 


-5.5094 


0.898655 


-61.669 


0.74580 


7 


44.34022 


1.000005 


12.38198 


0.948070 


-5.5096 


0.898654 


-61.666 


0.74570 


8 


44.34012 


0.999997 


12.38198 


0.948069 


-5.5095 


0.898657 


-61.669 


0.74581 


9 


44.34013 


0.999999 


12.38198 


0.948069 


-5.5096 


0.898656 


-61.670 


0.74580 


10 


44.34016 


1.000000 


12.38198 


0.948069 


-5.5095 


0.898657 


-61.669 


0.74582 


PVR 


44.34016 


1.000000 


12.38198 


0.948069 







4. We plan to investigate this problem in a forthcoming 
study. 

Further we shall include higher partial waves and em- 
ploy realistic NN forces to discuss physics. One expects 
that evidence for three-nucleon forces is more pronounced 
in the high energy region and the presented method is 
applicable there, now in the four-nucleon system. 



culations are performed on SX-5/128M8 (Research Cen- 
ter for Nuclear Physics), SX-5/6B (National Institute 
for Fusion Science), HP X4000 (Frontier Research Cen- 
ter for Computational Science, Tokyo University of Sci- 
ence) in Japan, and partly on Hitachi SR8000 (Leibnitz- 
Rechenzentrum fur die Miinchener Hochschulc) in Ger- 
many. 
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